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Abstract 

Decays of the D mesons to two pseudoscalars, to two vectors, and 
to pseudoscalar plus vector are discussed in the context of broken 
flavor SU(3). A few assumptions are used to reduce the number of 
parameters. Amplitudes are fit to the available data, and predictions 
of branching ratios for unmeasured modes are made. 
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Introduction 



Many data are available on the hadronic two-body decays of charmed mesons. 
Theoretical models that attempt to systematize the decay patterns have been 
available for many years. These models usually make dynamical assumptions 
in order to reduce the number of amplitudes that contribute to a particular 
decay. For example, the large N c approximation jl @, or the heavy-quark 
effective theory ||. It is not clear a priori how well such approximations 
should work and hence how seriously to take a conflict between a prediction 
and a measured value. Another approach is to assume that the matrix el- 
ements factorize 0]. This model is quite successful in describing observed 
modes, but again, it is difficult to know whether a discrepancy is due to 
an incorrect measurement of the failure of the assumption. A more general 
approach based on a diagrammatic classification ||, with different assump- 
tions also exists. In many cases attempts are made to obtain predictions of 
unmeasured modes from these models. 

577(3) is badly broken in these decays, so models based on exact sym- 
metry |J are not useful. An attempt at a complete parameterization of 
the 5*?7(3)-breaking has been conspicuously missing, due to the large num- 
ber of reduced matrix elements involved. We set out to remedy this omis- 
sion. This work gives a full parameterization of the decays of the D mesons 
into final states of two pseudoscalars (PP), two vectors (VV) and a pseu- 
doscalar plus a vector (PV), including S'[/(3)-breaking. The elements of 
this parameterization-the particle representations, the weak hamiltonian, the 
breaking operator, and the reduced matrix elements-are discussed in the fol- 
lowing sections. We make only very few assumptions to limit the number of 
parameters. We fit the parameters to the available data of two-body decays 
and predict many unmeasured modes. Because a few of the parameters are 
not constrained, we indicate which branching fractions are needed to predict 
the rest of certain classes of modes. We comment on the case of D s — > r]'p + , 
where the model is barely consistent with data. 
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1 Particle States in Flavor 577(3) 

In a model based on flavor SU(3), the particles are denoted by their £77(3) 
representation. The fundamental representation is the triplet (3) of quarks 
u, d, and s. The three D mesons {D°, D + , £>+} form an antitriplet (3) 
representation. The pseudoscalars {tt + , 7r°, 7r~, K + , K°, K~, K°, r] 8 } form 
an octet (8) representation, as do the vectors {p + , p°, p~, K* + , K*°, K*~ , 
K*°, us}. The rji and u)\ are each singlets. The physical r), r]', (ft, and u are 
linear combinations of them, with mixing angles -17.3° and 39° || for 
Tj-rj' and <p-uj respectively.^] 

2 The Weak Hamiltonian 

The decays of the D mesons are mediated by the weak hamiltonian. Ignoring 
QCD corrections, the hamiltonian in terms of the quark fields is 

#weak = % cos2 °c vrf(l - 7s)d s 7m (1 - 7s )c 

+ ^cos8 c sin6 c u^(l - 7 5 )s S7 M (1 - 7 5 )c 
- ^cos9 c sin6 c u^(l - j 5 )d d'y^l - 75 )c 
- ^sin 2 6 c M7 M (1 - 75)s d7 M (l - 7 5 )c. 

Note that the operators g create quarks and so transform as a triplet, while q 
transforms as the antitriplet. Using the Clebsch-Gordan coefficients for the 
expansion of the product 3 x 3 x 3, we can classify the operators according 

iK* denotes if* (892); r] 1 denotes 77' (958). 
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to irreducible representations of SU (3) as follows: 
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where (gg') denotes 97^(1 — 7s)g'- The numbers in parentheses are hyper- 
charge, total isospin, and third component of isospin of the particular mem- 
bers of the 577(3) representations. The weak hamiltonian can now be written 
in terms of the representations 3, 3', 6, and 15 as 
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Note that the 3 and 3' representations do not appear in the uncorrected 
J^weak [§]. Because the QCD corrections are multiplicative and do not mix 
the SU(3) representations, the 3 and 3' will also not appear in H wea ^(m c ). 

Since the decays of the D mesons occur at the scale of the c-quark mass, 
we must allow the QCD evolution of the various operators from the W- 
mass scale, where Equation ([!]) is valid, to the c-mass scale. The operators 
represented by the 15 are symmetric under quark interchange, and those rep- 
resented by the 6 are antisymmetric. The QCD renormalization of operators 
with these symmetry properties has been calculated flnj. We find that 
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in the regime where there are Nf flavor degrees of freedom. Taking into 
account the change in the number of active flavors as the 6-quark threshold 
is crossed, and using a(Mz) = 0.119, we obtain 

15 -> 0.81 15, 
6 -> 1.5 6. 
With Equation ([!]) as the boundary condition, we have 



weak 



sin 2 9 C 



+ cos 6 C sin 9, 



0.81 15(|,1,0) - 1.5 6(|,0,0) 
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(6) 



(7) 



Note that the QCD corrections do not introduce any new phases into the 
process. Unfortunately, until the values of the reduced matrix elements (dis- 
cussed below) are known, the coefficients in Equation (0) are of little use. 



3 Parameterization 

3.1 SU{3) Breaking 

For a complete parameterization of any process in flavor SU(3), we must 
include explicit breaking. Since we know that the source of flavor SU(3) 
breaking among the pions and kaons is the difference between the quark 
masses, we do this with an operator M which transforms as an 8. Although 
the quark mass difference is insufficient to explain the large SU(3) breaking 
that will be found, an octet is the simplest nontrivial operator that can be 
used. 

We can express M as 

M = a\ 3 + p\ 8 , (8) 

where A, are the usual Gell-Mann matrices. The term in a represents break- 
ing of the isospin SU(2) subgroup. This breaking, proportional to the differ- 
ence between up and down quark masses, is expected to be very small and we 
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neglect it in the following. The constant (3 can be absorbed into the reduced 
matrix elements. Hence M can be reduced to 

M = A 8 . (9) 
3.2 Reduced Matrix Elements 

Now consider the most general parameterization of the decays in the context 
of the flavor 577(3) symmetry. For each possible contraction of the represen- 
tations into an SU(3) singlet there must be one parameter, i.e., one reduced 
matrix element. The representations involved are those in Section 1: D (3), 
H (6 + 15), and two of P and V (each 1 or 8). In addition, we must include 
all possible ways of involving the symmetry-breaking parameter M. We as- 
sume that the breaking is linear in M. Each reduced matrix element is, in 
principle, complex. We have chosen to contract D with H, then contract 
the products (PP, PV, VV) (and then possibly with M), and finally contract 
the two parts into the singlet. Our labels for the reduced matrix elements 
reflects this. For example, the matrix element denoted (DH 15 ) 8 ((PP) 1 M) 8 
is obtained by contracting D and the 15 component of H into an octet, con- 
tracting PP into a singlet which combines with M to become another octet, 
and contracting the two resulting octets into the singlet. 

Unfortunately, the above parameterization involves far more parameters 
than there exist data. Therefore we make two important assumptions. First, 
we assume that we can separate the spin and flavor dynamics of the pro- 
cesses, i.e., that the relative strengths of the reduced matrix elements are 
the same in the PP, PV, and VV cases. This implies that only forty-eight 
reduced SU(3) matrix elements are needed. They are labeled with S and 
O for the singlet and octet representations, rather than with PP, PV, or 
VV. In order to distinguish the spin states we introduce two parameters, 
called "PV" and "VV." Second, we assume that the phase of each reduced 
matrix element is given solely by the representation of the product parti- 
cles (before M is included). Bose symmetry for PP and VV and an ap- 
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propriate phase rotation of the particle fields reduces the list of phases to 
(T^Oi, (7^)1, (P Vl ) 8 , (PwOb, (Vr/Og, (VwOb, (PP)x, (PP) 27 , 

(PV)i, (PV>, (PV)io, (PV)r , (PV) 27 , (W)i, and (VV) 27 . One should 
note that we cannot determine the relative phases between PP, PV, and VV. 

The amplitude for each decay mode can be expressed as a sum over the 
reduced matrix elements with the appropriate Clebsch-Gordan coefficients: 

A{D j - Xi) = J2C ijk R k Si. (10) 

k 

Here R k are the reduced S77(3) matrix elements and Si are the parameters 
that we call PP = 1, PV, and VV. The SU(3) Clebsch-Gordan factors C l 
were calculated by computer. Many of the routines used are described in 



3.3 Linear Combinations of Reduced Matrix Elements 

There are 45 measured values for the two-body decay modes and an addi- 
tional 13 modes where upper limits exist .| It would appear that there are 
still more parameters than data, and therefore the model lacks predictability. 
However, there are only forty linearly independent combinations of the 577(3) 
reduced matrix elements that contribute to the possible decay modes of the 
D mesons. With the assumption of the last section concerning the phases of 
the reduced matrix elements, the linear combinations fall into these classes: 

§The data are from the Particle Data Group |8), together with |L3| for the mode D + — > 
K*°n+. 
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We write them each as a sum over the reduced matrix elements, viz., 

% 

and normalize them for convenience by setting 

Ec£ = i- (12) 

i 

Now Equation fllCf ) is replaced by 

A(D J ^X l )=Y,C'; jn L^S l . (13) 

n 

The It"* repl ace the reduced matrix elements in our parameterization of 
the amplitudes. The forty linearly independent combinations contain matrix 
elements including those that involve the breaking operator M. It is not 
possible to divide the linear combinations into a set that contains only matrix 
elements without M and a set containing only matrix elements with M. Of 
the forty combinations, three are not constrained by the available data. We 
call them L^ 2 \ and L^ 3 \ They are discussed below. 

The replacement of the set of reduced matrix elements by the set of 
linear combinations that contribute to the possible decay modes reduces the 
number of parameters by eight. The total number is now fifty three. These 
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parameters are fit to the data; the individual reduced matrix elements are 
no longer considered. 

The unconstrained combination contributes to the modes D° — > rjrj, 
r/i]', r/'r]', r)<f), tjuj, r/'cj), rfui, rfxfi, 4>lj, and ojuj. Because these modes are un- 
observed, the phases of (771^1)1, and {uiiUJi)\ are also unconstrained. The 
remaining unconstrained linear combinations are and L^K They con- 
tribute to the above modes, and also to modes of the types D° — > r/K and 
D s — > r]K + . By "type" we mean a class of modes that contain mesons of 
the same flavors and charges. Thus the type D s — > i]K + contains the modes 
D s -> r]K + , r]'K + , r]K* + , r]'K* + , 0K + , uK + , <pK* + , uK* + , and no others. 
With the exception of the limit on the branching ratio for D s —>■ 4>K + , there 
are no data for these modes. We still have some freedom in the definition of 
and that allows modes of the type D° — > r]K° to depend on only one 
of them (choose L^). This will allow us to estimate one of their branching 
fractions and thereby make some predictions of the other modes of this type. 



4 Data and Fitting Thereof 

The data used to determine the parameters are listed in Tables 1-5. These are 
the modes for which there exist either experimental values or experimental 
limits. In the VV modes, S and D waves are possible. Data exist from E691 
14] for the modes D° K*°p° and D + — > K*°p + . These are consistent with 
the S- and D-waves both having significant amplitudes and are inconsistent 
with either being zero. The ratios of S- and D-wave amplitudes from these 
two modes are taken as additional data, and the overall ratio of S- to D- 
wave amplitudes for the VV modes is allowed to vary in the fit. Its value is 
determined by the two modes mentioned above, and depends very little on 
the other data. 

For each mode we remove the phase space and Cabibbo factors and reduce 
the branching ratio to a decay amplitude in arbitrary units. Because the 
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vector particles have substantial widths, the phase space for modes involving 
a vector is integrated over the relativistic Breit-Wigner for that resonance. 
The effect of this is important for those modes where the sum of the particle 
masses is within a few widths of the D-mass. The modes D° — > <j)K*°, <j)K*°, 
and D + — > <pK* + would be forbidden if the widths were set to zero. Each 
amplitude is now expressed as a sum of Clebsch-Gordan coefficients times 
the parameters that represent the reduced matrix elements, and finally as a 
sum over the linearly independent combinations of reduced matrix elements. 

The parameters were fit to the data amplitudes with MINUIT, release 93.11 
[ E| . The total \ 2 was found to be 30.9 for seven degrees of freedom, indi- 
cating that the overall fit was poor. However, more than half of the x 2 arose 
from only one mode. The mode in question is D s — ► i]'p + . The experimental 
value for the branching ratio D s — > r]'p + cannot be accommodated in our 
scheme. It is measured [15| to be larger than that for D s — ► r]p + , an a priori 
surprising result. We note that the angular distribution of the decay pions 
is barely consistent with that expected. A confirmation of this experimental 
value would be very significant as all other models [|J also predict a ratio of 
B(D S — > i]'p + )/B(D s — > f]p + ) of less than one. 

We decided to reject the experimental value for the branching fraction of 
D s — ► r]'p + . The result is a better fit, from which the branching ratios are 
reported in the tables. The total \ 2 is now H-6 for six degrees of freedom. 



5 Predictions 

5.1 Predictions from the Fit Parameters 

From the fit values of the parameters the branching ratios of decay modes 
were calculated. In Table 1 are presented the modes for which there exist 
experimental values. Our calculated branching ratios are consistent with 
the data, with the exception of D + — > K*°K* + and D s — * t]'p + . For the 
former the fit prefers a branching ratio that is three standard deviations 
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below the reported experimental value. The latter was removed before the 
fit (see Section 4) because its experimental value was questioned. For this 
mode we predict a branching ratio of (1.5 well below the reported 

experimental value ||15|| . Tables 2-4 contain modes for which there is no 
experimental information or for which there is an experimental limit. We 
have attempted to predict the branching ratio of each mode from the fit. 
However, in some cases the uncertainties are so large that we are able only to 
provide (90% confidence level) limits on the branching ratios. Notice that in 
all cases in which there are experimental limits, our predicted branching ratio 
or predicted limit is in the allowed region. We are unable to say anything 
about the mode D s — > p°ir + , because the uncertainty on its prediction is 
greater than the experimental limit. 

There are two modes, D s — > 7r + 7r° and D s — > p + p°, which are forbidden 
in a model without isospin breaking. They are predicted to be identically 
zero. The modes that are kinematically forbidden are D° — > rfrf, rj'<p, and 
(fxf). The modes involving the linear combinations lP~\ L^ 2 \ and are 
discussed below. Any PP, PV, or VV mode not appearing in the tables is 
higher order in the weak coupling Gp. 



5.2 Unconstrained Linear Combinations 

There remain three linearly independent combinations of the reduced matrix 
elements that are not constrained by the data. The combination con- 
tributes only to modes of the type D° — > 7777. contributes to the types 
D° — > 7777 and D° — > r]K°. contributes to these modes, and to modes of 
the type D s — > rjK + . 

The first unconstrained linear combination LW contributes only to ampli- 
tudes involving (SS)\. These amplitudes, it is worth noting, are due entirely 
to SU (3) breaking. However, when we include the phases, we must make four 
estimates in order to obtain two predictions of modes of the type D° 7777. 
This would be an unproductive endeavor, and so we forego it. 
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In order to predict the modes of the types D° — > r\ K° and D s — > i]K + , 
we need two new inputs. In order to show the variability of the resulting 
predictions, we try three different sets of inputs. Scheme A is motivated by 
the recent CLEO measurement of the doubly-suppressed mode D° — > n~K + 



17], in which this mode is found to have a branching ratio of about three 



times that of the corresponding unsuppressed mode, D° — > n + K . For this 
scheme, the two inputs are 



B(D° -> 7]K°) = 3 tan 4 6 C B(D° -> rjK°) 
B(D S -> 77^+) = 3 tan 2 8 C B(D S -> ^tt+) 



(14) 



The linear combinations L^ 2 ^ and Z/ 3 ) are then constrained and the remaining 
branching ratios in the column for scheme A in Table 5 are found. The 
predictions for scheme B are based on the following estimates: 

B(D° -> r]K°) = 3 tan 4 6 C B(D° -> r]K°), 
B(D S ^4>K+) = 3tan 2 fl c B{D S -> 0tt+). ^ ^ 

A third scheme (C) is considered also. It is based on these estimates: 

B(D° -> = 3 tan 4 fl c £?(_D° -> 0^°), 

B{D s ^<pK + ) = I tan 2 #c B(D S — > </>7r + ). ' ' ° 

The resulting predictions are again in Table 5. The spread in these values 
provides an indication of the expected ranges for these quantities. 

One should note that arbitrary choices of the above modes may fail to give 
an acceptable fit, given the constraints from measured modes. For example, 
an apparently reasonable choice would have been 



B(D° -> r]K°) = tan 4 6 C B(D° -> rjK ) 
B(D S -> r]K + ) = tan 2 6 C B(D S -> ?77r+) 



(17) 



A consistent fit cannot be obtained to implement this. The parameters 
and could not be given values to accommodate B(D° — > r)K°) < 
0.0052%. 
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5.3 Modes Involving Axial Vectors 

There are a few modes involving axial vectors that have been observed or for 
which there are experimental limits. However, those that involve if(1270) 
and -?T(i400) are mixtures with the l + ~ octet, which we can call B since 
it includes the b\ (1235). Therefore, in order to include these modes in our 
framework, we require two new parameters, "PA" and "PB." In addition, 
we must also accommodate the mixing between fu x \ an d /i(8) to become 
/i(1285) and /i(1510), as well as the new phases that are introduced. There 
are too few experimental observations of the PA and PB modes to make this 
endeavor fruitful. For that reason, they are not included here. 



6 Comments on Models 

It is clear from the data alone that significant SU(3) breaking is necessary 
in any successful model of D decays. For example, B(D° — > K + K~) = 
B(D° — > 7r + 7r~) in exact SU(3), yet they are in reality quite different. Models 
based on exact SU(3) @, |§, ]18[ (or even on nonet symmetry [|19j) are thus 



not admitted by the data. 

Models of D decays based on heavy-quark effective theory (e.g., ||) have 
as yet not developed to the point at which individual nonleptonic decays can 
be calculated. The question of whether HQET is applicable to the c quark 
is still unsettled. The HQET is based on an expansion in the parameter 

^-0.2 (18) 

m c 

and assumes that it is small. Certainly this would be a good assumption 
in the case of the b quark, but perhaps not so here. Until we are able to 
calculate branching fractions in HQET, we must reserve judgement on its 
applicability to the D mesons. 

Diagrammatical methods to the problem of D decays present us with a 
complementary approach to the one adopted in this work. The parameters 
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in the SU(3) framework represent sums of diagrams in the diagrammatical 
approach. A very general diagrammatical calculation of branching fractions 
appears in || . Two shortcomings of their work lie in final-state interactions 
and in the inclusion of £77(3) breaking. The phases of the final-state inter- 
actions are added to the model, and are external to its central theme, and 
therefore appear as an ad-hoc mechanism to force a fit. SU(3) breaking is 
added to the calculation as an additive correction to the diagrams in which 
it is believed to be important. However, there is also hidden breaking in the 
addition of phases in the final-state interactions. The result is a model in 
which the size and source of ££7(3) breaking is not easily discerned. It is 
difficult to draw any conclusions from the application of such a model. 

The factorization method is a special case of the diagrammatical ap- 
proach. In it certain diagrams are considered unimportant (i.e., the anni- 
hilation diagrams). However, ||] find that these diagrams must be again 
included, as well as final-state transitions and intermediate resonances. The 
result is an eclectic model with little elegance. We are unable, because of 
the ad-hoc features, to comment on the reliability and predictability of this 
model. 

A description of nonleptonic D decays in a large-iV c (number of colors) 
expansion || is an elegant one with few parameters. In it, the source of 
££/(3)-breaking is introduced my including nearby resonances. It is also a 
subset of the diagrammatical approach and neglects some diagrams based on 
their suppression by 1/N C . One may argue that these diagrams are larger 
than thought, and cite the fit of || as evidence of this. Nevertheless, [|| 
obtain excellent agreement with the data, with the exception of some modes 
involving rj and r]'. In this model, SU(3) breaking is introduced only through 
the inclusion of resonances in one class of diagram. They obtain, in agreement 
with our work, large breaking. 
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Conclusions 



There now exist enough data to constrain all but three combinations of the 
reduced matrix elements of the broken SU(3) model of the decays of D 
mesons with the two assumptions discussed in Section 3.2. We have used 
these data to do so. Using the experimental information on 57 modes we are 
able to predict branching ratios or upper limits for an additional 53 modes. 
Only two measured modes are not easily accommodated in the fit. The 
measurement of a few additional modes involving r], r] f , 0, uj would enable 
another dozen or so modes to be predicted. 
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Table 1: Modes with positive experimental values. Branching ratios from 
data and from the fit are given. 



Mode 
~W 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 





Data BR 


Fit BR 


K- 7T + 





0401 


± 





0014 





0400 


± 


0.0014 


K- K + 





00454 


± 





00029 





00453 


± 


0.00030 


K° 7T° 





0205 


± 





0026 





0208 


± 


0.0022 


K° K° 





0011 


± 





0004 





00103 


± 


0.00043 


7T~ 7T + 





00159 


± 





00012 





00159 


± 


0.00012 


7T- K + 





00031 


± 





00014 





00031 


+0 

— 


00018 
00014 


7T° 7T° 





00088 


± 





00023 





00087 


± 


0.00025 


ri K° 





0068 


± 





0011 





0069 


± 


0.0011 


rf K° 





0166 


± 





0029 





0168 


± 


0.0028 


K*~ p + 





059 


± 





024 





063 


± 


0.016 


K*° p° 





016 


± 





004 





0164 


± 


0.0038 







0029 


± 





0015 





0029 


+0 
-0 


0019 
0014 


u> K*° 





011 


± 





005 





0099 


± 


0.0044 







0019 


± 





0005 





00192 


± 


0.00045 


K~ p + 





104 


± 





013 





102 


± 


0.013 


K~ K* + 





0034 


± 





0008 





00323 


± 


0.00080 


K° p° 





0110 


± 





0018 





0110 


± 


0.0017 


K*~ 7T + 





049 


± 





006 





049 


± 


0.0058 


K*~ K+ 





0018 


± 





0010 





00209 


± 


0.00087 


K*° 7T° 





030 


± 





004 





0301 


± 


0.0039 


<f)K° 





0083 


± 





0012 





0081 


± 


0.0012 







020 


± 





004 





0195 


± 


0.0043 







019 


± 





005 





0204 


± 


0.0049 
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Mode 



Data BR 



Fit BR 



D + 
D + 
D + 
D + 
D+ 
D + 
D + 
D + 
D + 
D + 
D + 
D + 



K° 7T + 

K° K + 

7T° 7T + 
7] 7T + 

K*° p+ 
K*° K*- 
K° p+ 

7T+ K*° 



K*° 

'*0 



7T 



K*° K + 

(j) 7T + 
fST+ 



0.0274 

0.0078 

0.0025 

0.0075 

0.021 

0.026 

0.066 

0.00046 

0.022 

0.0051 

0.0067 

0.00039 



± 0.0029 
± 0.0017 
± 0.0007 
± 0.0025 
± 0.014 
± 0.011 
± 0.025 
± 0.00015 
± 0.004 
± 0.0010 
± 0.0008 
± 0.00020 



0.0262 

0.0086 

0.00257 

0.0068 

0.0398 

0.0090 

0.071 

0.00046 

0.0217 

0.00463 

0.00674 

0.00039 



± 0.0028 
± 0.0016 
± 0.00067 
± 0.0021 
± 0.0092 

+0.0054 
-0.0041 

± 0.018 
± 0.00014 
± 0.0041 
± 0.00097 
± 0.00078 

+0.00027 
-0.00020 



D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
D, 



K° K + 

K° K* + 
K*° K+ 

(j) TT + 



V P 



+ 



0.035 


± 


0.007 


0.0319 


± 


0.0059 


0.019 


± 


0.004 


0.0204 


± 


0.0039 


0.047 


± 


0.014 


0.054 


± 


0.012 


0.056 


± 


0.021 


0.055 


± 


0.018 


0.065 


± 


0.017 


0.056 


± 


0.014 


0.042 


± 


0.010 


0.043 


± 


0.011 


0.033 


± 


0.005 


0.0328 


± 


0.0053 


0.035 


± 


0.004 


0.0349 


± 


0.0040 


0.100 


± 


0.022 


0.100 


± 


0.019 
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Table 2: D° modes with predicted branching ratios. Experimental limits are 
given when available. All limits are at 90% confidence. 



Mode 



Data BR 



Predicted BR 



Predicted limit 



D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 
D° 



o K o 



K*° 
p° K*° 

7T° p° 
P~ 7T + 

VP° 

K *- K *+ 

p° K° 
rf K*° 
p- K* + 
p- K+ 

7T° 
TT" K*+ 
7T° K*° 

K° K*° 
ujtt 
r}' p° 

LOp° 
p- p + 

K*° K° 
rf 7T° 
p° p° 



< 0.0011 



< 0.0008 



< 0.0015 



0.00017 

0.00108 

0.00038 

0.00014 

0.093 

0.094 

0.0060 

0.025 

0.0024 

0.0024 

0.00018 

0.00022 

0.0020 

0.024 

0.0019 

0.0020 



+0.00011 
-0.00008 

+0.00073 
-0.00054 

+0.00031 
-0.00022 

+0.00018 
-0.00011 
+0.133 
-0.075 

+0.136 

-0.076 

+0.0092 

-0.0050 

+0.041 

-0.021 

+0.0041 

-0.0021 

+0.0041 

-0.0021 

+0.00032 

-0.00016 

+0.00038 

-0.00020 

+0.0035 

-0.0018 

+0.049 
-0.022 

+0.0037 
-0.0018 
+0.0041 
-0.0019 



< 0.092 

< 0.0092 

< 0.0091 

< 0.00070 

< 0.00085 

< 0.0078 

< 0.0080 

< 0.0087 

< 0.00052 

< 0.086 

< 0.011 

< 0.084 

< 0.015 

< 0.00061 

< 0.057 

< 0.0065 
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Table 3: D + modes with predicted branching ratios. Experimental limits are 
given when available. All limits are at 90% confidence. 



Mode 



Data BR 



Predicted BR 



Predicted limit 



D+ 
D + 
D+ 
D + 
D + 
D + 
D+ 
D + 
D + 
D + 
D + 
D + 
D + 
D + 
D + 
D + 
D + 
D + 
D + 
D + 
D + 
D + 
D + 
D + 



7] K + 

7T+ K° 
7T° K + 
7T° p+ 

K*+ 
p+ K*° 
p° K*+ 
p+ K° 

00 K* + 

UJ 7T + 
7T° K* + 

UJ p + 

K° K*+ 

V P + 
rf K + 
p° K + 

1]' 7T + 

uj K + 

1 p + 

p° 7T + 

V K* + 
T] K* + 



< 0.007 



< 0.012 



< 0.009 

< 0.015 

< 0.015 

< 0.0014 



0.0066 

0.0032 

0.017 

0.0086 

0.0034 

0.00031 

0.025 

0.0095 

0.0087 

0.0022 

0.0024 

0.0103 

0.0026 

0.0012 

0.0012 

0.0016 

0.0018 

0.00094 



± 0.0023 

+0.0030 

-0.0020 

+0.018 

-0.011 

+0.0089 

-0.0057 

+0.0036 

-0.0023 

+0.00035 

-0.00022 

+0.031 

-0.018 

+0.0118 

-0.0071 

+0.0119 

-0.0068 

+0.0031 

-0.0018 

+0.0036 

-0.0020 

+0.0162 

-0.0087 

+0.0049 
-0.0023 
+0.0026 
-0.0011 
+0.0022 
-0.0011 
+0.0041 
-0.0015 
+0.0042 
-0.0017 

+0.00237 
-0.00092 



< 0.011 

< 0.0054 

< 0.0048 

< 0.0082 

< 0.0086 

< 0.0048 

< 0.0074 

< 0.0012 

< 0.00071 

< 0.00091 

< 0.000082 

< 0.0022 
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Table 4: D s modes with predicted branching ratios. Experimental limits are 
given when available. All limits are at 90% confidence. 



Mode 



Data BR 



Predicted BR 



Predicted limit 



D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
D s 
Ds 



7T° K + 
7T + K*° 

i p + 

7T° K* + 

p° K*+ 
p+ K° 
p° K+ 
uj p + 

7T+ K° 

K° K* + 
K° K + 

p + K*° 
K*° K+ 

UJ 7T + 
7T° P + 
p° 7T + 
7T+ 7T° 



12.0 ± 3.0+ 



< 0.007 



< 0.017 

< 0.0028 



0.0059 

0.038 

0.015 

0.077 

0.0126 

0.031 

0.049 

0.012 



= 
= 



+0.0048 
-0.0034 

+0.047 

-0.028 

+0.019 

-0.011 

+0.096 

-0.058 

+0.0164 

-0.0096 

+0.043 

-0.024 

+0.071 
-0.040 

+0.030 
-0.012 



< 0.061 

< 0.0015 

< 0.00039 

< 0.00046 

< 0.00057 

< 0.0080 

< 0.00025 

< 0.0090 

< 0.064 

not significant 



tSee text. 
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Table 5: Modes based on estimates. The only available experimental limit is 
shown. Values marked with * are inputs. 



Mode 




Data BR 


Fit BR 

(scheme A) 


Fit BR 

(scheme B) 


Fit BR 

(scheme C) 


D° -> 






0.000054* 


0.000054* 


0.00035 


D° -> 


V K° 




0.00046 


0.00046 


0.00085 


D° -> 


(j) K*° 




0.000019 


0.000019 


0.000016 


D° -> 


uo K*° 




0.00027 


0.00027 


0.0012 


D° -> 


(j) K° 




0.00054 


0.00054 


0.000066* 


D° -> 


uK° 




0.000096 


0.000096 


0.0014 


D° -> 


T] K*° 




0.00048 


0.00048 


0.00094 


D° -> 


v > K *o 




0.0000083 


0.0000083 


0.0000024 


D s - 


T] K+ 




0.0027* 


0.00041 


0.0031 


,0, - 


i K + 




0.017 


0.052 


0.015 


,0, -> 


K*+ 




0.011 


0.024 


0.0095 


,0, -> 


u K*+ 




0.0057 


0.028 


0.0046 


£> s - 


(j) K+ 


< 0.0025 


0.00051 


0.0033* 


0.00037* 


£> s - 


u K+ 




0.0064 


0.019 


0.0055 


,0, - 


T] K* + 




0.00083 


0.00015 


0.00094 


D s - 






0.00090 


0.0028 


0.00077 
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